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We investigate the dynamics of quantum particles in a ratchet potential subject to an ac force 
field. We develop a perturbative approach for weak ratchet potentials and force fields. Within 
this approach, we obtain an analytic description of dc current rectification and current reversals. 
Transport characteristics for various limiting cases - such as the classical limit, limit of high or low 
frequencies, and/or high temperatures - are derived explicitly. To gain insight into the intricate 
dependence of the rectified current on the relevant parameters, we identify characteristic scales and 
obtain the response of the ratchet system in terms of scaling functions. We pay a special attention 
to inertial effects and show that they are often relevant, for example, at high temperatures. We find 
that the high temperature decay of the rectified current follows an algebraic law with a non-trivial 
exponent, j oc T~"^^. 



PACS numbers: 05.40.Jc, 05.60.Gg, 73.23.Ad, 85.25.Dq 



I. INTRODUCTION 

Ratchets have attracted a considerable recent inter- 
est because of their paradigmatic role as microscopic 
transport devices (for review articles, see e.g. Refs. 
|l|,||j3|). Applications range from microacale electron- 
ics - includingpthc photogalvanic Ci^ectjj transport in 
quantum 



junctionsi 




and antidot arraydJ, - over Josephs^ 
and vortex mattert2l to cell biology. 
At the same time, ratchets are of fundamental theoretical 
interest since they represent one of the simplest nonequi- 
librium systems. 

The analysis of ratchet systems reaches back quite 
some time before Feynman drew the attention of a 
wide audience to such systems in his lectures where 
he discussed the possibility to employ ratchets as heat 
enginesO Subsequently, researches in ratchets have been 
progressing steadily, in parallel in different scientific com- 
munities, until the explosive outbpist in theoretical and 
experimental interest in the 1990al3 had occurred. 

In this paper, we report on the analytic progress in the 
study of the so called tilting ratchets, where the combina- 
tion of an asymmetric static potential with an unbiased 
ac force and a coupling to a heat bath leads to current rec- 
tification. Past theoretical studies pf Ais ratchet type fo- 
cused on the classical massless casetjO and have revealed 
the current reversal phenomenon, i.e., the possibility that 
the direction of the rectified current reverses its direction 
when model parameters such as the frequency or ampli- 
tude of the ac current are changed. The inclusion of a 
finite mass of the particles has shown -that it may give 
rise even to multiple current reversals.EJ Further exten- 
sions accounted for the quantum nature of particles and 
of the bath. Quantum fluctuations were fpttBd-tp provide 
an additional source of current reversals .t3E2lEil 

In essence, the direction and amplitude of the current 
turned out to be very sensitive to the various system 



parameters. While this dependence makes ratchets valu- 
able for applications - such as devices that can separate 
particles of different species - it still lacks a satisfying the- 
oretical understanding. Analytic approaches can give in- 
sight into this problem. However, even the single-particle 
problem is already so complex that analytic approaches 
can be adyanced only in limiting cases,, such as the adia- 
batic limitllj or the deterministic limit .E3l3 



In our paper we develop a perturbative approach valid 
for weak ratchet potentials and weak driving forces, 
which covers a wide range of practical applications. 
Within this perturbative approach we are able to capture 
all prominent phenomena including multiple current re- 
versals. This approach provides a unified framework for 
deriving and understanding the dependence of the recti- 
fied current on the particle mass, temperature, friction 
coefficient, and frequency of the driving force. We pay 
particular attention to the role of inertial effects and show 
that they lead to a substantial current enhancement even 
in the high- temperature limit. 



In Sec. II we specify the model and establish a path- 
integral formulation as analytic framework. The pertur- 
bative scheme is developed in Sec. [IT In Sec. [V we 



briefly demonstrate that the linear mobility can be con- 
veniently obtained from this approach and that results 
for special cases known in the literature are reproduced. 
However, ratchet effects can be obtained only in nonlin- 
ear response. The leading nonlinear mobility is calcu- 
lated and evaluated for various limiting cases in Sec. 0. 
We conclude with a discussion of our approach and re- 
sults in Sec. VI. Technical details of our calculations are 



presented in appendices. 
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II. MODEL 

We consider a quantum particle of the mass to in a sta- 
tionary ratchet potential U{x). In addition, we impose 
an ac driving force F{t) which is chosen to be unbiased, 
i.e., to vanish-jupon time averaging. Following Caldeira 
and LeggettJ^ we couple the particle linearly to a bath 
of harmonic oscillators at temperature T. This bath si- 
multaneously provides friction and a fluctuating force for 
the particle. For simplicity, we assume a linear spectral 
distribution of these oscillators, giving rise to Ohmic dis- 
sipation. In the classical limit, the particle coordinate 
x{t) follows the equation of motion 



mi{t) = -U'{x{t)) + F{t) - r^x{t) -t- ^{t) 



(1) 



with a friction coefficient 77 and a Gaussian thermal noise 
f (i) obeying 



(2) 



To account for the quantum nature of the particle and 
of the bath, we follow the analiajisi-of quantum Brownian 
motion by Fisher and ZwergeicJCJ who studied the case 
of a sinusoidal potential and a dc driving force. 

The rectified particle velocity V can be determined 
from the average particle coordinate X{t) via 



V 



X{t) 



hm -Xit), 
t— ►00 t 



dx X P{t, a;). 



(3a) 
(3b) 



where P(i, x) is the probability distribution for the par- 
ticle position at time t. This distribution is related to 
the reduced density matrix operator p{t) (after the bath 
degrees of freedom are traced out) by 



P{t,x)^{x\p{t)\^) 



(4) 



[We use the Dirac notation, where {x'^\p{t)\x^) is the 
density matrix in position representation.] 

The dynamics of this density matrix is most conve- 
niently treated in . the Feynman- Vernon path integral 



representation.LjIfll The time evolution of the density ma- 
trix from some initial time ti to a final time is given by 



{xj\p{ti)\Xf ) JJ dx+dx. J{tf,x+,Xf ;ti,x^,x. ) 

x{x+\piU)\xr) (5a) 

with the kernel 

J{t{,xl ,xY;ti,x^ ,x^) = J Dx Dy (5b) 

being a double path integral over a|U.|trajectories x{t) and 
y{t) with the boundary conditionsE3 

xiti,i) = ^K+f + yiti.f) = T{xt,f - x^'^i)- (6) 



The path integral involves the effective action 

S = Sq -\- Si, 

So ^ I ff dtdt' v{t)K{t - t')v{t') 



(7a) 



+i J dt y{t)[mi{t) + r]x{t)], (7b) 

Si ^ tj dt[Y,Ys^ix{t) + sy{t))-y{t)F{t)],{7c) 

s 

with all time integrals running from ti to tf. For no- 
tational convenience, the usual contribution U{x'^{t)) — 
U{x^{t)) is written as a sum over the spin-like variable 
s — ±fi./2 which, however, does not have the meaning of 
a physical spin. 

The effective action (^ includes already the aver- 
age over the bath degrees of freedom. This average 
leads to an integral kernel Kit) which reads in Fourier 
representatiorO (we set fee = 1) 



K{ijj) = rjhjj coth 



(8) 



In the classical limit, K{ijj) — 277T reproduces the corre- 
lator (^). For T — K{uj) — rih\uj\ represents a kernel 
that is highly nonlocal in time representation. 

The model has a large number of parameters: the par- 
ticle mass m and the friction coefficient 77, then h and 
T as a measure of the strength of quantum and thermal 
fluctuations. Further parameters are implicit in U{x) 
and F[t). The potential can be represented by a Fourier 



U{x) = Y,U,e^'^ 



(9) 



with amplitudes Uq for wave vectors q. For periodic po- 
tentials with a period a, the wave vectors are 



27r 

q = n — 
a 



with integer n. 

In analogy to U, the ac drive is represented as 



Fit) 



(10) 



(11) 



with Fq = since the force is assumed to be unbiased on 
time average. For a periodic drive with period tp, the 
frequencies uj are integer multiples of the basic frequency 
2Ti/tF- Although we assume here periodicity of U and F, 
a generalization to random U and F is straightforward 
and will be discussed in the end of this paper. 



III. PERTURBATIVE APPROACH 

The definition (||) of the velocity has the drawback 
that one has to calculate Xit) as the expectation value 
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of the final position in an ensemble of forward-backward 
paths of a finite length if — t;. In order to avoid techni- 
cal complications related to boundary effects, we relate 
the average velocity to an expectation value at an inter- 
mediate time t which can be kept fixed while the limit 
ti —oo and cxd is been taken. 
Consider the "partition sum" 



Z = dx 



dx^dx-^ J{tf,Xf,Xf;ti,x^ ,x^ ) 



x{x+\p{ti)\x. 



(12) 



of all forward-backward paths between U and tf. It is 
normalized to Z = 1 since it is the trace of the density 
matrix at t{. We define 



V{t) ^ {x{t)) 



(13) 



as expectation value in the ensemble Z of fluctuating 
paths. In this definition, we can take the limit U — oo 
and tf ^ oo right away. In the absence of a nonequilib- 
rium driving force and due to the presence of dissipation, 
V{t) would vanish after an initial relaxation for every 
possible initial density matrix p(ti). In the presence of 
the driving force and in the limit U — oo, V{t) will be 
determined uniquely by F{t) and independently of the 
initial state. 

AlthO|Mt|We strictly follow the definition of Fisher and 
Zwergerr^ri in the path integral formulation of the prob- 
lem, we differ in the definition of the average velocity. We 
argue in appendix ^ that, in the long time limit, the time 
average of the velocity V{t) coincides with the earlier def- 
inition (|^) in combination with (^ and (^. We find the 
expectation value V{t) to be a convenient quantity for 
the subsequent perturbative evaluation. 



A. Perturbative expansion 

To make analytic progress, we consider F as small and 
calculate the nonlinear dynamic response of the velocity 
to the driving force. 



V{t) 



dt'ni{t~t')F{t') 

dt'dt"n2{t -t',t- t")F{t')F{t") 



1 

2! . 
-OiF"^). 



(14) 



The mobilities pm can be expressed conveniently as ex- 
pectation values in the path ensemble using the partition 
sum as generating functional. 



pi{t-t') = 



SV{t) 



5F{t') 



F=0 



tM2{t-t' ,t-t") = 



= (iW*y(i')>lF=o' 

5^V{t) 



F=0 



SF{t')6F{t") 
{x{t)zy{t')zy{t")) 



F=0 • 



(15a) 



(15b) 



The generalization to higher order mobilities is straight- 
forward. The expectation values now refer to the equilib- 
rium system in the absence of the driving force. 
After Fourier transformation, Eq. (|l^ ) reads 

Vuj = pi{uj)F^ + ^ fi2{i^',uj")Fuj'F^»Suj,uj'+u," 



-OiF"^). 



(16) 



The rectified current is given by the time-average (zero 
frequency component) of the velocity, 



Since the driving force is unbiased, Fq — 0, current rec- 
tification cannot be obtained in linear response. Rather, 
ratchet effects require frequency mixing which is present 
only in nonlinear response. For weak F the leading 
ratchet effect will be determined by fi2- 
If the driving force has the symmetry 



F{t)^-F{t-to) 



(18) 



for some time to (for example, if F is monochromatic), 
the rectified velocity will be invariant under the transfor- 
mation F{t) F{t - to) = -F{t). Then the contribu- 
tions to the rectified current from all mobilities /i™ with 
odd m must vanish. 

Although the calculation of these mobilities is already 
much simpler than a closed calculation of V{t), it still 
cannot be performed analytically for general potentials. 
Therefore, we employ a second expansion in U , utilizing 
the weakness of the potential. 

The mobilities, which, according to Eqs. (E5|), are the 
equilibrium expectation values jim = {Om)\F^Q of ob- 
servables Om = i{t)iy{t') ■ ■ -iyi^^^), will be calculat 
perturbatively in the potential using the expansionEJI 




Tl = 



^ 2is 



We thus can write 



.(19) 



(20) 



n=0 



with 



h, E nil/*-/*- 

n 

exp {% ^ \x{t,) + s,y{t,y^ )^ (21) 



,(") - 



where the average (• ■ ■ ) o is governed by the "free" action 
So defined by Eq. (7b). Since 5*0 is Gaussian, the av- 



erages can be performed straightforwardly using Wick's 
theorem. 
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B. Free theory 

For these averages it is important to know the corre- 
lations of the free theory. In Fourier representation, one 
easily finds 

{x{lj')x{uj"))o = C{lo')S{uj' + (22a) 
{x{u')iy{Lu"))a = GiLu')S{Lo' +Lo"), (22b) 
{yiu;')yiu"))o = 0, (22c) 



We conclude this subsection by pointing out some key 
features of the response and displacement function. For 
h = 0, G and W are related through the fluctuation- 
dissipation relation 



(28) 



In the quantum case with T — 0, W{t) diverges for all t 
in the limit m —^ 0. 



with the response and correlation functions 
G(^) ^ 



irjuj + muj"^ — 0+ ' 
K{u:) 



(23a) 
(23b) 



To calculate the nonlinear mobilities, one has to use the 
retarded response function which is in time representation 



G(t) = (a;(t)iy(0))o = -[1 
with a relaxation rate defined by 



m 



e-^']Q{t) (24) 



(25) 



Causality of response 

G(t) = for t < 



(26) 



is reflected by the Heaviside step function 8(t) in (24). 
Note that {x'^{t)) = oo since the free system is translation 
invariant and the particle spreads diffusively (subdiffu- 
sively for T = 0) over the entire space. Therefore, C(t) 
is not a well defined quantity. Instead, the displacement 
function 



W{t) 



1 



(Ko-x(o)no 



duj 
2^ 



[1 - cos{ijjt)]C{uj) 



(27a) 
(27b) 



C. Characteristic scales 

Before we move on to a further evaluation of the path 
integral, we pause for a moment to fix the relevant time, 
length, and energy scales of our problem. From the re- 
sponse function of our problem we can identify the typical 
relaxation time 



_ 1 _ 



7 



■q 



(29) 



Rewriting the displacement correlation function W{t) — 
PW(t/tr:e\) in tcrms of the dimensionless function W of 
the dimensionless argument t/t^ci, we identify from Eqs. 
( ^ ) the diffusion lengths i for the thermal and the quan- 
tum case. 



^2 



Tm 



^2 



h 
v' 



The de-Broglie wave length 



niT 



(30) 



(31) 



is a related further characteristic scale for the particle in 
the absence of dissipation. 

Alternatively to Eqs. (^0[), we can associate with ther- 
mal and quantum fluctuations a characteristic energies 

E = rj^i'^/m, 



Eth = T, Equ = ^7- 



(32) 



captures all information about the particle 

(sub)diffusion. This quantity will play a central 

role in perturbation theory. Unfortunately, it can be 
calculated explicitly only in limiting cases: 



'77 

W{t) ^ln7|t| for r = 

TTT] 

For semiquantitative purposes. 



1] for h = 0, (27c) 
(27d) 



W{t) « ^[7|t|-|-e-^l*l 



l] + ^m + (7in (27e) 



The potential and driving force - which act as probes to 
the free particle - define the space period a, time period 
tp, and amplitudes 



U = 2\Ua 



T = 2|F„ 



(33) 



defined by the lowest harmonic modes q and uj. In the 
case of random U or F, the periods would be replaced 
by a correlation length or time and the amplitudes by 
variances. 

In terms of these scales, a necessary requirement for 
the validity of the perturbative approach is that exter- 
nal probes must be weak in comparison to the internal 
fluctuations, i.e.. 



is a good interpolation over the whole parameter range. 



U, aT <C max(£'th, i?qu) 



(34) 
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These scales will also determine the location of the phe- 
nomena under consideration, as we will discuss later on. 
However, as we recall by calculating the linear response 
mobility, th condition (|34| ) is not sufficient for the validity 
of perturbative results. 

In the subsequent calculations it is convenient to use 
dimensionless quantities. It is natural to choose as time 
scale, or 7 as frequency scale. The generic length scale 
is the potential period a. The ratios of ^^j^ and i'^^ to 

provide a natural measure for the strength of thermal 
and quantum fluctuations. Hence we define 



i = "ft, (2) 



1 



-, qj=aqj, 



as dimensionless quantities. 



D. Mobilities 



In Eq. (^), the mobilities /im^ are determined by 
expectation values which can be calculated conveniently 
from a generating functional. We define 



exp I 



(^i / dT[p{T)x{T) + (7{T)y{T)] 



(36) 



as functional of auxiliary fields p(r) and (t{t). r is a 
real time variable and we will be distinguishing it from 
t only for bookkeeping purposes. The averages deter- 



mining the mobilities /i- 
functional derivatives 



(n) 



can be then represented as 



where one has to identify 



Cr(T) = ^gjSj(5(T-tj), 



(38a) 
(38b) 



after performing the functional derivatives. Using the re- 
sults of the previous subsection, the generating functional 
can be expressed as 



Z = exp / / dridr; 



-p(ti)C(ti - T2)p{t2) 



-ip(ri)G(ri - T2)<t{t2) 



(39) 



As mentioned previously, C{t) is divergent. This implies 
that Z = if / dTp{T) = J2j 1j 0- Therefore, Z can 
be nonvanishing only if the "momentum conservation" 
1j = is satisfied. In this case one may rewrite 



Z — exp 



// 



dTldT2 



rl 



p{ti)W{ti - T2)p{t2) 



-ip{Ti)G{Ti - T2)(t{t2) 



(40) 



exp {i ^ qj[x{tj) + Sjy(ij)]) 



d ,5(™+i) 



' dt5p[t)6a{t')---5a{ti'^')y 



(37) 



Hereby, we introduce the abbreviation J^qj = ^Y^qjfi- 
The subsequent calculations of the mobilities are based 
on this generat ing functional. For later convenience, we 
combine Eqs. (^l|), (^7|), and (Eol) to our master formula 



. d 1 



si.qi.,--- .s„,qr, j=l ' '' 



dtr, 



(5/9(i)(5cr(t')---^o-(i^™^) 



exp (^JJ dTidT2 ^^p{ti)W{ti ~ T2)p{t2) + ip{Ti)G{Ti - T2)(j{t2] 



(41) 



Thereby, the substitution ( p8| ) has to be made after all 
functional derivatives are taken. Momentum conserva- 
tion implies that all mobilities vanish for no ?i = 1. For 
n = 2 and even m the mobilities vanish since the con- 
tributions to the sum in the right hand side of Eq. (|4ll ) 
are odd in {q}. We have noted already above, that no 
contribution to the rectified current can arise from fim 
with odd m and arbitrary n if the driving force obeys 
the symmetry (fsh. In this case, up to fifth order in F 



and U, the only contribution comes from 1I2 



(3) 



After we have determined the generating functional Z 
for the mobilities in the previous section, we now turn to 
the evaluation of the lowest order mobilities of interest. 
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IV. LINEAR MOBILITY m 

Although we do not expect ratchet effects from linear 
response, it is instructive to calculate in order to 
verify that the present calculation of the mobility repro- 
duces that results of Fisher and ZwergeilHSEa for static 
F and sinusoidal U (i.e., for this purpose we include the 
amplitude -Fb in our consideration). 



A. Leading orders 



To zeroth order in U, it is obvious that 



To first order, 



(42) 



(43) 



since the momentum conservation mentioned above can- 
not be satisfied (strictly speaking, it is satisfied for the 
mode q — which, however, does not enter the dynam- 
ics). 

To second order, a straightforward calculation (see ap- 
pendix ^ leads to 

^,['\u;) = iu;G\u;)Y,q'\U,\'AB^^lju.) (44) 



with 



AB^'ljcj) ^ B'_^l^{u; = 0) - B'_^l^{cj), (45a) 
= ^sin[^9?G(0]e-''?^W<5,,+,„o.(45b) 
In the last expression, the sine has a nonunique sign for 

(46) 



= r(2) 



(2) 



n<-q G(oo) = — 



reflecting quantum interferences of particle trajectories. 

We briefly discuss interesting limiting cases of ^\ for 
which we will also examine the ratchet effect later on. 



B. Classical limit 



with 7(-, •) the incomplete gamma function (to be dis- 
tinguished from the parameter 7). We introduced the 
dimensionless frequency 



_ .2 2 Tmq^ 
'-'q = 4h9 = 



(49) 



related to the thermal diffusion time over a distanc e 1/1 7 
via Witj-ci/jyo) = q~^. The insertion of expression (48b) 
into Eq. (Bj) yields an explicit analytic expression for the 
classical linear response mobility at finite frequencies. 



1 



1 



7] (1 — iu})'^ 

q 



l{.Vq - i'^-, Vq) (50) 
for the special 



which reproduces Eq. (4.11) of Ref. 
case of sinusoidal U and oj — Q. 

In the massless (overdamped) limit m ^ 0, where B 
can be easily Fourier transformed, this simplifies to 



^^2^1 ^ T 

with a dimensionless scaling function 



(2) 



/ 9 



j2 i^^i' 



IZ 



(51) 



(52) 



Thus, for m = 0, the corrections to mobility decay pro- 
portional to T~2 at high temperatures. 



Adiabatic limit 



The adiabatic limit w — > simplifies the calculation of 

,(2) 



AB_^^^{uj), resulting in 

,(2) 



Mi'(^-O) = -^T.\Uq\ 



X / dtt e-?'^(*)sin[-g2G(t)](53) 



In a classical limit, h 0, both the overdar nned p 
underdamped cases are understood fairly well.l23E2lE3'E3 
In the present perturbative approach, the fluctuation- 
dissipation theorem (ESh allows for the simplification 



1 d 



(47) 



In this case, the Fourier transformation can be performed 
analytically. 



(48a) 



which agrees with the linear response limiting case (4.18) 
of Ref. ||. 

At r = 0, the particle shows a remarkalJeJocalization 
transition due to the dissipative coupling.LJi2J For strong 
coupling, the particle is localized in an arbitrarily weak 
potential, whereas it remains mobile for weak damping. 
This transition is reflected by the divergence of the mo- 

(2) 

bility correction n\ {uj = 0) due to a divergence of the 
time integral at large t. From the logarithmic asymp- 
totics (27d) of W{t) one can identify the location of the 
transition at a = 1 with 



1 



27rn 2TTh 



(54) 
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Note that for a < 1 the inequahty (|4^) is fulfilled for all 
wave vectors q, i.e. quantum interference effects suppress 

(2) 

the contribution to ^\ . 

In the strong damping regime, the divergence of 
^\ {u! = 0) signals the break-down of perturbation the- 
ory. Thus, at r = the condition a ^ 1 should be added 
to the condition (p^). This condition may be regarded 
also as a condition for the period of the potential (with 
localization for < 2tt£'^^). 

In the perturbatively accessible regime of a < 1, Fisher 
and ZwergerO pointed out the interesting fact that the 
mobility is a nonmonotoneous function of temperature. 
At zero temperature, the particle has its free mobility. 
Weak thermal fluctuations {T < T*) first reduce mobility 
(thermally resisted quantum tunneling), whereas strong 
thermal fluctuations increase mobility back to its free 
value (thermally assisted hopping). The crossover occurs 
for a <C 1 at the temperature 



(55) 



at which the de-Broglie wave length is comparable to the 
potential period, A ~ a. 

Before we continue to enter new territory, we wish to 
conclude this subsection by stressing that our approach 
successfully reproduces previous linear response results 
for Lj — Q and already provides additional insight into 
the frequency dependence. 



V. NONLINEAR MOBILITY ^2 

The generating functional formalism presented above 
provides an efficient tool to calculate also higher order 



mobilities. Here, we focus on the lowest order of ^2 in U 
contributing to current rectification. 



A. Leading orders 



To zeroth order. 



M^"^ {t~t\t~ t") = {x{t)iy{t')iv{t"))o = (56) 



vanishes since S'o is invariant under the reflection 
{x, y\ — > {—a;, To first order, 



(57) 



vanishes again because momentum conservation cannot 
be satisfied. To second order, 



^l^2\t - 1' - 1") ^ Q 



(58) 



according to the general statements following Eq. (^l|). 

The general third order contribution {t — t' ,t ~ t'') 
is given by expression (C6) calculated in appendix 0. 
Since this expression is somewhat clumsy and since we 
are interested only in ratchet effects, we can restrict our 
considerations to 



^ Uq.Uq^Ug^qi 



X {q,q2 [2i?f;, (0, 0) - Sj (-C., 0) - b5 {lo, 0)] 



'{q} 

-qi<Z3[2i3j5(0,0) 

j(3) 



(3) 



-Kq2<?3[2i3|;;(o,o) - i?j;;(o, -u) - i3|;;(o, w)]} 



(3) 



(59) 



with 



B 



{^}(*l -t2,t2- h) 



^ sin[^giGi2g2] ^ sin[^ {qi Gisqs + 92623^3)] 

X exp(gi 14^1252 + 92^^2393 + 9lW^1393)'5gi+?2+93,O023- 



(60) 



Note that is an implicit function of {q} and invariant 

f 31 

under {q} — + —{q}- Consequently, changes 
sign under a reflection Uq — > U-q which implies that the 



rectified velocity ( p7| ) vanishes for even potentials as it 
should. Examining the contribution from a set of wave 
vectors {q} and its reflected set —{<;}, one can recognize 
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that {—1^,1^) is real and that it depends only on 



with a reduced scaling function 



^qi Uq^ 



(61) 



-2 E t^S^i+^s) 



(In order to obtain simple analytic expressions below, we 
find it most convenient to use momentum conservation 
to eliminate the sum over 52-) In Eq. ([59|), the sum over 
the momenta can be restricted to qj 7^ (in addition to 
momentum conservation) since B^^^ vanishes otherwise 

(later on, these restrictions are referred to by Yl')- Phys- 
ically, it is clear that a constant shift of the potential 
cannot enter the dynamics of the particle. 

Eq. ( p9| ) is our main result in general form. Further 
analytic progress is hampered by the absence of an ana- 
lytic expression for W{t). Nevertheless, further analytic 
progress is possible in various limiting cases. 

Using the dimensionless quantities defined in Eqs. ( |35| ) 
we may reexpress Eq. (p9) as 



fif(h,f,^) (62) 



with jl being a dimensionless function of dimensionless 
arguments. For a monochromatic driving force, the rec- 
tified velocity is 



fi-^^ {h, f, u) 



to leading order according to Eq. (|6|). 

B. Limit h and m 



(63) 







We start with the examination of the classical limit. To 
further simplify the analysis and to perform a comparison 
of our perturbative results with previous approaches, we 
consider the overdamped limit with m = 0. In this case, 
the memory kernel B) { simplifies considerably to 



qiq2ql 1 1 
\q} 2 '■ '■ — 



with the characteristic frequencies 
Vq = Tq^/rj. 



(64) 



(65) 



This frequency corresponds to the time a. classical 
particle needs to diffuse over a distance q~^ ■ Insertion of 
Eq. (p|) into Eq. leads to 



(66) 



(67) 



2gi _ gig3(gigi - 

qf + z^ {qf + z^)iql + z^) 



(68) 



The primed sum is restricted to momenta satisfying mo- 
mentum conservation and qj ^ 0. It is interesting to 

realize that the function /t2ci(-^) uniquely determined 
by the shape of the potential. If current reversals ex- 
ist in the limit under consideration, they correspond to 
oscillatory behavior of this function. 

The scaling function /i2'^ci(^) becomes simple in addi- 
tional limiting cases. In deriving these limits from Eq. 
( |68| ) , one has to make use of momentum conservation and 
of permutations of momentum labels. For 2 — > cxd 



(3) -3;;3- 
9l92 93 



91 «2 93 



dx [U"'{x)YU'{x). (69a) 



The terms of order cancel each other. Thus, we easily 
retrieved the result obtained previously in Ref. ^ 
In the opposite limit z ^ 0, 



2 y m- 



2 



with a potential integral 



91 

dx U^{x)6T{x) (69b) 



T(x) 



dy Uiy), 



ST{x) = T{x) ~ i / dyT 
a Jo 



(70a) 
(70b) 



In deriving Eq. (69b) we have assumed Uq = 0, otherwise 
additional subtraction terms should be added. 

This scaling behavior ( |69| ) implies the asymptotic be- 
haviors of the rectified velocity 



Vq oc T u! for uj — > CXD, 
Vo oc lu"T-* for T ^ 00. 



(71a) 
(71b) 



The apparent divergence of Vq for T — > is an artifact 
of leaving the range of validity of our perturbative ap- 
proach. Analogously, the apparent divergence of Vo for 
77 is due to the assumption of overdamped dynam- 
ics. Nevertheless, these divergences may be interpreted 
as indications that ratchet effects are particularly strong 
at low T and in the underdamped case. This situation 
will be examined later on. 
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T 



(3x^ ^ r-, FIG. 2: Vo(T) for oj = 0.01 (bold line), oj = 1 (long dashes), 

FIG. 1: Scaling function ii^^^,{z) for the potential (|7g) shown ^ ^ 4 (^j^^^.^ dashes), u = 5.5 (dotted), and = 7 (dash- 
in the mset. dotted) for comparison with Fig. la in Ref. |l^ The vertical 

line represents the vicinity of the current reversal for uj = 1. 



Before we move on to other limiting cases, we show 
that the current reversal phenomenon is captured by our 
perturbative approach. We find it also instructive to per- 
form an explicit quantitative comparison of our results 
with the exact results in Re f. | l7| . For this comparison, 
we evaluate Eqs. (^ with (|67D for 



U{x) = 



sin(27r-) + - sin(47r-) 
a 4 a 



(72) 



cf. the inset of Fig. The corresponding scaling func- 
tion ( |68| ) is shown in Fig. |^. Since it has one zero, we 
expect one current reversal. 

We explicitly compare our perturbative result for the 
rectified velocity as a function of temperature - displayed 
in Fig. H - with the exact solution displayed in Fig. 
la of Ref. |l^. Thereby, length, energy and time scales 
are fixed by the choices a — \, U — and rj — 1. 
The monochromatic driving force is F{t) — T sm{Ljjt) 
with amplitude T = 0.5. From the shape of the scaling 
function it is clear that we find a current reversal with 
varying temperature for every lo > and also a current 
reversal with varying frequency for every finite T. The 
quantitative agreement is good for T ^ U, where the 
perturbation theory in U is justified. 

The classical limit is not restricted to single current 
reversals. It is likely that an arbitrary number of current 
reversals can be obtained by suitably tailored potential. 
We have found, for example, that it is sufficient to add 
one more harmonic to obtain a second current reversal. 
Specifically, the potential 



0.02 



0.01 



<=!. 




100 200 300 400 



500 



FIG. 3: Scaling function fi^lilz) for the potential ( [t^ shown 
in the inset. 



C. Limit T (X> for m > 

It is interesting to examine the high-temperature limit 
since there are significant differences between the cases 
m = and m > 0. For T — > 00, the exponential factor in 
Eq. (|60| ) strongly suppresses and thus fj,'^\ In this 
limit, we find the asymptotic behavior of the mobility 
(see appendix ^ 



1 



' UJ 



-(3) 
2 ^2,hT' 



(74a) 



U{x) = -U 



sin(27r— ) H — sin(47r— ) H — sin(67r— ) 
a 4 a 4 a 



(73) 



or 



leads to the scaling function shown in Fig. ^ with two 
zeroes, i.e., two current reversals. 



2 \ 17/6 
77a 7 ^ 



.(3) 

1 + Lo'^li- 



:(74b) 
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with 



-(3) 



9i>0,?3>0 



r>(3) ~2 



(91 + 93)^93 



1/3 



31/3^2 



- (3) 

The constant is uniquely determined by the shape 

of the potential. In the high temperature limit, the rec- 
tified current is much stronger for massive particles than 
for massless cases, since ^'"2 ^ ~ j^-n/e for m > whereas 
^(3) _ rj.-4 fQj. ^ = 0, cf. Eq. ( |71b| ). This observa- 
tion is consistent with the mass dependence in Eq. (74), 

fj,^^ ~ nJ I^T^^'^ 1^ for TO ^ 0, which signals that in the 

(3) 

limit TO ^ 0, [I2 should decay with a higher power of 
temperature. Thus, inertial terms are crucial at high 
temperatures even for large friction where the relaxation 
of the particle in the minima of the potential is over- 
damped. 




log,n [1/(113-^7)] 

FIG. 4: Double logarithmic plot of the rectified velocity 
[the dimensionless quantity fi^^]^ given in Eq. ([tg])] versus 
temperature in the adiabatic limit for the potential in 
the underdamped case Q = | (bold line). The dashed line is 
a guide to the eye representing the behavior oc of the 

high-temperature limit. 



D. Limit to 



For large frequencies, Eq. (Efl) simplifies to 



1 



W2(l+w2) 



(75a) 



or 



M2^^(-w,w) 



r;(?7a2w)2;i2(^2+^2)A^2,M 



{h,f) (75b) 



with 



Ull\a'qi{ql + qiq3 + ql} 



(75c) 



[For the discussion of this hmit, i?|^|(0, 0) = _B|^j(tj = 

f3^ 

0,w = 0).] ^2 iif is a function of the potential shape 
and of to parameters measuring the strength of quan- 
tum and thermal fluctuations. In the special case h = 
m — Q we found in Sec. 



dence i?j'^i(0, 0) oc 92/91 which led to a cancellation in 



VB a momentum depen- 



the sum over momenta in expression (75c). Using the 
fluctuation-dissipation relation (|2^) , one can easily show 

0. Thus, this 



that _Bj^j(0, 0) is independent of to for h - 



0, i.e., A2^hf(0'^) 



0. 



cancellation persists as long as h 

For TO > this implies a decay ^'^\—uj,oj) oc uj^^. How- 
ever, such a cancellation can no longer be expected for 
h > 0. In this case, one finds again ^i^^"^ {—i^, 1^) oc uj~^ at 
large frequencies. 



E. Limit u ^ 

A further limit of interest is the adiabatic limit for 
the quantum particle. This limit was also studied in the 
pastll3'E2ltll and has revealed that additional current re- 
versals may arise from the competition of quantum and 
thermal fluctuations. 

For — > 0, Eq. (|5^ ) reduces to 



(3) 
{9} 



dt'dt" 



•^mii' -qd"fBf^^{i\i") (76a) 



with 



- fe2n(3) 



(76b) 



We have calculated i^^2\d numerically for h ^ ^ [i.e. 
a = rjo? / (2'kK) = ^ corresponding to the delocalized 
case, cf. Eq. (p4[)] as a function of temperature (cf. Fig. 
^ for the potential The two poles of the double- 

logarithmic plot in Fig. H represent current reversals. At 
high temperatures, the relation ji^2^ oc y-i^/e jg recov- 
ered (dashed line). At zero temperature, a finite current 
is generated by quantum fluctuations. 



VI. CONCLUSIONS 

We have developed a perturbative approach for quan- 
tum ratchets, which captures current rectification and 
reversals of the current direction. Our main results are 
the analytic expression ( p9[ ) for the leading nonlinear mo- 
bility and its evaluation for various limiting cases. In 
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particular, the high-temperature Hmit for massive par- 
ticles has revealed the relevance of inertial terms even 
for strong damping. Since the rectified current decays 
like Vq oc T"^ for massless particles whereas it decays 
like Vo oc T-^^^ for massive particles, inertial effect can 
lead to a substantial enhancement of ratchet effects. On 
the other hand, in the high-frequency limit, the quan- 
tum nature of the particle is important. While Vq cx 
for massive classical particles, quantum fluctuations also 
enhance the rectified currant, leading to Vq esc uj~'^. 

While our perturbative approach is limited to weak po- 
tentials and driving forces, it has the advantage that it 
can be easily generalized to higher dimensions. There- 
fore, applications for example to asymmetric antidot 
arraysQ become possible. Furthermore, a generalization 
to random ratchet potentials is obvious. Thereby one 
could describe the case of asymmetric potential wells with 
random positions. This generalization can be achieved if 
one allows for continuous wave vectors q of the potential 
and simply replaces Uq-^ Uq^ Uq^ by its average in the non- 
linear mobility (^9|). An extension of this perturbative 
approach from single quantum particles to electron gases 
is under current investigation by the authors. 
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APPENDIX A: VELOCITY 

In this appendix we show that the average velocity can 
be calculated from the definition (13). It is this definition 
in which we. deviate from the approach of Fisher and 
Zwerger.blEl 

The original definition (|3^), in a more explicit form, 
is based on the average distance the particle travels in a 
large time interval. 



V 



Xjtt) - XjU) 

ti - u 



(Al) 



in the limit t{ oo and — > — oo. Thereby the position 
expectation value ( |3b[ ) and the time evolution ( pa[ ) of the 
density matrix. On the other hand, the time average of 



J 



V{t), Eq. (|13[), can be written as 

X{tf)-X{U) 



V{t) 



ti-ti 



(A2) 



with 



X{t) = J dxxP{t,x), (A3a) 

P{t,x) = {S{x-x{t))). (A3b) 
A priori, P(t,x), which is an expectation value in an 
ensemble of paths of length tf— ii, is different from P{x, t) , 
Eq. (^), which is an expectation value in an ensemble of 
paths of length t—ti. However, the definitions coincide for 
t = ti and also for t = ti. In the first case the definitions 
coincide, in the second case because one can integrate 
out the paths (the integral corresponds to Z, Eq. (|T^), 
the integral is most easily performed for a diagonal initial 
density matrix). Thus, 



X{ti) - X{ti) - X{tt) - Xiti) 



(A4) 



and V = V{t). If there is a well defined expectation value 
V{t) = {x{t)) = ■^X{t) for ti — oo and tf oo, it must 
coincide with V{t) since boundary effects from times near 
ti and tf should become negligible in this limit. 



APPENDIX B: CALCULATION OF fif'> 

Here we present intermediate steps of the calculation 
leading to Eq. (|4|). In a first step, we need to evaluate 

S^^Swf ^ {^Git-t') + JJdr,dMWit-n)pin) 

+ iG{t - Tl)a(Tl)]p(T2)^G(T2 - t')}z (Bl) 
where Eqs. (^sj) have to be inserted for n = 2. Thereby, 

2; ~ ggiWl2<32+9liGi2S2<32+92iG2lsigijJ ^ ^ (B2) 

where we abbreviate Wm = W{tk — ti) etc. and we used 
W{Q) = 0(0) = 0. Note that in the last exponential, 
6*12 or 6*21 vanishes for all ii, t2 because of causality. 
Inserting Eq. (Bl) into Eqs. (36) and ( [2l]), only the 
last of the three terms coming from Eq. ( p31| ) survives 
summation over Sj in Eq. (pT|). One obtains 



t') = 



^ 4siS2 

<3l92SlS2 



dtidt2G{t - ti)siqi[qiG{ti - t') + q2G{t2 - t')]Z. 



(B3) 
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In these remaining terms, the summation over Sj yields 

^ (t-f) = / / dtidhG{t - h)qi [qiG{h - t') + q2G{h - (*1 " ^2), (B48 



9192 



^(2)^(^1-^2) = ^^siZ-|sin(Jg2gi2)cxp(-g2T4^i2)<5,,+,„o. (B4b) 



Fourier transforming this expression leads to Eq. (HJ). 



APPENDIX C: CALCULATION OF 



(3) 



Following the same route as for /i^^-* , we first calculate 



/)3 

Z ^ I drp(r) [iG{t - t')iG{T - t") + iG{t - t")iG(r - t')]-2^ 



(5p(t)(5cr(t')^o-(i") 



+ j dT[W{t - t)p{t) + iG{t - r)cr(r)] 

X / dr' p{t')iG{t' - t') f dr" p{T')iG{T" - t")Z. (CI) 



Eq. (H) leads to 

2 = g9l 1^1292+91 W'i3l33+q2Vl^2393 + [92iG21+93*G3l]9lSl+[9ljGi2+g32G32]92S2 + [gi2Gi3+g2i^ ^ (C2) 

with (5{q} = <5qi+ij2+g3,o- Considering the right-hand side of Eq. ( pl| ) as a sum of four contributions, the first three 
disappear after summation over the Sj. For example, if ti < t2 < t^, Z is independent of S3 and the summation over 
S3 leads to a cancellation. The remaining fourth contribution to Eq. ( pl| ) reads explicitly [623 = 0(^2 — ^3)] 

.^^^(t-t',t-.") . - E kkklll "^'^^ e-G(.-t,)sr.r 

9l,92,93,Sl:S2,S3 ^ .i J J J 

x{q,G{h - t')qiG{h - t") + q2G{h - t')q2G{t2 - t") + q:,G{h - t')q:iG{h - t") 
+qiq2[G{ti - t')G{t2 - t") + G(<i - t")G{t2 - t')] 

+gig3[G(ii - t')G{h - t") + G{h - t")G{h - t')] 

+9293[G(i2 - t')G{h - t") + G(<2 - t")G{h - t')]}Z (C3) 

where we used permutation symmetries among indices j which allow us to restrict the time integrals to t2 > t^. Then, 
summation over {s} leads to 

B^.,^(t,-t2,t2^t^ ^ 623 E ^^^^^^^^ (C4) 

S1S2S3 

= e23;^sin[JqiGi2g2]sin[^(giGi3g3 + 92G23g3)]e«^'^-«^+«^'^=^«^+'^^^-'^5{,}. (C5) 
Note that ^jg} 7^ only for ti > t2 > t^. In terms of sj^j = Bqflj^q^ we obtain 

Hi^\t-t\t-t") ^ t E Uq,Uq,Uq, f f f dhdt2dt3 G{t - tl)qi 
91,92,93 ■^■^■^ 
x{q,Gih - t')qMh - t") + 92G(t2 - t')q2G{t2 - t") + q^G{t^ - t')q^G{h - t") 
+9i92[G(ii - t')Git2 - t") + G(<i - i")G(t2 - t')] 
+9i93[G(ii - t')G(t3 - t") + G(<i - i")G(t3 - t')] 

+q2qAG{t2 - t')G{h - t") + G(<2 - i")G(t3 - - ^2,^2 - ^3) (C6a) 
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or, 



after Fourier transformation, 

lj!i\uo',uo") = [uo' +uj")G{uj' + u")G{u:')G{uo") ^ U^.U.^U^.qi 

qiq2 93 



Ratchet effects are related to 



9192 93 

X {q^ijJJ) (0, 0) + q^^sf^j + 0) + qg^sf^j (c.' + c.'', c.' + c.'') 
q2qMl\ iuj' + uj'\ UJ') + B[l\iio' + Lu", lu")]}. 



Uj" ~ —LO' = LO, 



(C6b) 



= LO, for which Eq. ( p9| ) foUows after usage of momentum conservation. 



APPENDIX D: DETAILS FOR T ^ oo 

Although straightforward, the calculation for the high-temperature limit requires some care. For this calculation, 
it is convenient to rewrite Eq. (Is^) as 



1 



7] rj'^LU^ + m^w^ 



J2 U,,Uq,Ug,qi dh2dt23F[l\{u;,h2,t23)exp[-E[l\{h2,t23)] (Dl) 



wfth tjk = tj - tk, trs = ti2 + t23, 



i^|^j(w,ti2,t23) = 2{gi(3'2[l - cos(wii2)] + qiqsll - cos(w<i3)] + (72'3'3[1 - cos(wi23)]} 

2 h 2 h 

x-sm[-qiG{ti2)q2]TSin[-{qiG{ti3)q3 + 920(^23)93)], 



h '2 



'h '2 



(D2a) 



and 



Ef^yiti2,t23) = -biM^(ii2)g2 + g2T^(i23)'?3 + giW(ii3)'?3] = ^(bia;(ti)+(Z2a:(i2)+'?3a;(<3)]')o>0. (D2b) 

I 



With increasing T, the rectified current shrinks since 
(3) 

E^gj increases proportional to temperature. The dom- 
inant contributions come from small ti2 and small ^23- 



We proceed with an expansion of E^^'y and F^^y 



(3) 



to ex- 



tract the leading orders for large T. 

In the high-temperature limit, one can neglect the 
quantum contribution to W and expands for small times 
(using q2 = —qi — qs because of momentum conservation) 



{9} 



Tm / (g| + g|)^r2 , 1 2 91+93 r3 
< ;2 + ''-^ 



772 [ 



93 



We introduced dimensionless times t^. and t_ via 

91 



7^12 
7*23 



-t. 



93 



93 



(D3) 

(D4a) 
(D4b) 



To extract the asymptotics for T ^ 00, one has to 
distinguish the contributions for qi/q^ > and qi/qs < 
(remember that one needs to consider only qi qs)- 



Because of causality, the time integrals cover only the 
quadrant with ti2 > and ^23 > in the (^12,^23) plane. 
This quadrant corresponds to ranges 

i+ > and - —i+ < f_ < —i+ for — > 0,(D5a) 
93 qi 93 



f_ > and — £_ < <+ < 
93 



for ^ < 0(D5b) 
91 93 



(D6) 



The integrals are transformed via 

dti2dt23 = ^^-^-l^di+di^. 

For 91/93 < 0, it is sufficient to retain the quadratic 
term in Eq. ( p3| ) since it implies that £_ ^ T^^l"^ . Then 
also t+ ~ according to Eq. (Dq), i.e., t\2 ^ t23 ^ 



T ^1'^ . Since is quartic in small times, 
_,2 i„. _L 

n^T - ' 95 93 

-0(T-5/2), 



^{9} ~ „2^,2 



(D7) 



the resulting contributions to [1^2^ will be of order T^^ . 
These terms can be neglected in comparison to terms of 
order T'^"^/^ which come from 91/93 > 0. 
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For qi/q-i > it is not sufRcient to retain the quadratic 
term in Eq. (D3) since the integral over would diverge. 



(3) 

Thus, one has to include cubic orders in E^^j which imply 



that i+ ^ T"i/3. Consequently, the higher order terms leadmg order 
not explicitly written in Eq. (D3) can be neglected. Since 

The 



[ignoring condition (D5)] since the quadratic term in (D3) 
provides a cutoff that dominates over the condition (D5) 
(the errors decay exponentially in T). Therefore, the 
iH+ - r-ii/6 will not result in a 



(3) 



(3) 
{■/} 



T ^ we now expect /ij'^^ 

A3) 



contribution to ^2^' of order T since it is odd in i_ 
Performing the time integrals for the remaining terms. 



proper expansion of Fj^-J up to order T now yields 



1 



E 



F, 



(3) 
{9} 



^ Q?, 



^2 
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- ^{11 + l3)P-ri} + 0(T-"/6). (D8) 

Thereby it is sufficient to retain even orders in i_ because 
the integral over t_ can be extended to all real values 



— ill 



qi,q3;qi/q3>o 
q3){qi+qlr I dU 



93 



91 



(D9) 



yields Eq. (|7|). 
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